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Abstract 

A density-functional theory is established for inhomogeneous su- 
perfluids at finite temperature, subject to time-dependent external 
fields in isothermal conditions. After outlining parallelisms between 
a neutral superfluid and a charged superconductor, Hohenberg-Kohn- 
Sham-type theorems are proved for gauge-invariant densities and a 
set of Bogolubov-Popov equations including exchange and correlation 
is set up. Earlier results applying in the linear response regime are 
recovered. 
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Experiments on confined Bose-condensed gases have revealed a rich va- 
riety of dynamical behaviours. These include elementary excitations of low- 
lying shape deformation modes propagation of sound waves in both the 
condensate and the thermal cloud 0, Josephson-like oscillations in double 
condensates 0, creation of vortices 0, phase dynamics from various atom- 
laser configurations and Bloch oscillations of a condensate in an optical 
lattice g]. 

While mean-field theories suffice in most cases to describe the observed 
behaviours, condensates can now be created where effects beyond mean-field 
can be explored, by tuning almost at will the scattering length and hence 
the condensate self-interaction energy or else by approaching the criti- 
cal temperature for Bose-Einstein condensation. Along these lines a Time- 
Dependent Density- Functional Theory (TD-DFT) for superfiuids is a suitable 
framework to treat their dynamics with inclusion of exchange and correlation. 

For a many-electron system in the normal state the foundations of the 
theory come from a set of theorems by Runge and Gross P, |^, which have 
been extended to superconductors at zero temperature by Wacker, Kiimmel 
and Gross (WKG) ||lOl- In applications to normal electron systems in the 
linear response regime, the limitations to low-frequency phenomena have 
been conceptually overcome by Vignale and Kohn [0]. Their current-density 
formulation of TD-DFT embodies plasmon dispersion and damping as well as 
transverse-current fluctuations, allows a unified treatment of the damping of 
collective excitations from the Landau and mode-coupling mechanisms and 



yields microscopic generalized-hydrodynamic equations 12 . 

A similar scheme has already been developed for the dynamic linear 
response of superfiuids |]13|, extending to inhomogeneous systems and to 



finite-frequency phenomena Landau's hydrodynamic equations in the two- 
fluid model. The present Letter concerns the underlying foundations of TD- 
DFT for superfiuids. The proof of the relevant Hohenberg-Kohn-Sham-type 
theorems parallels the WKG derivation for superconductors and one of our re- 
sults is a "dictionary" which translates vector and scalar potentials. Maxwell 
equations and the like from a charged to a neutral fluid characterized by 
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spontaneous symmetry breaking. We also allow for finite temperature in 
isothermal conditions. This preludes to a specific choice of the reference 
system for the TD-DFT mapping, which is described by a set of Bogolubov- 
Popov equations including the non-condensate density. 

Introductory material. The dynamics of a system of interacting spinless 
bosons confined in a static potential and evolving from an initial equilib- 
rium state at time to is driven by the Hamiltonian 

H{t) = fj^it) + Vv{t) + S,{t) + W{t) . (1) 

The system is subject to an external vector potential A and to scalar fields 
V and r]. We have Ta = -{l/2m) J drip^{r,t) [hV - iA{r,t)f '^{r,t), Vy = 
J drV{r,t)i>^r,t)^lj{r,t) and S^, = J dr \r]{r,t)ijj^r,t) + r]*{r,t)ip{r,t)] [0]. 

W{t) is the interaction term, given in terms of the field operators by W{t) = 
(1/2) / dr^dr2iP^r^, t)tP^r2, t)w{r^, r2)V'(r2, t)V'(ri, t). 

A couples to the total current density and V to the total particle den- 
sity. After writing the ensemble average of the field operator as {ip{r,t)) = 



nc(r, t) exp[i(p{r, t)], we see that the symmetry-breaking source field rj drives 
both the density of condensate nc(r, t) and its phase ip{T,t); the latter de- 
termines the irrotational part v<j(r,)f:) of the velocity field through \s{r,t) = 
{h/m)'Vip{r,t). In the linear regime S*^ we can write S*^ = / dr[X{r,t) ■ 
5vs(r, t)+a{r, t)6nc{r, t)] in terms of the condensate-density operator 6nc{r, t) = 
2Re[('?/'(r, to))5^"l'(r, t)] and of the irrotational- flow operator (5vs(r, t) = {h/m)'V60 
with 60{r,t) = — Im[5-?/'t(r, t)/ (^/'"''(r, to))]. The (real) fields a and A are then 
related to rj by a{r,t) = [nc{r,tQ)]^^Re[{ilj'^{r,to))ri{r,t)] and V ■ A(r, t) = 
-2mIm[(7/'t(r,to))77(r,t)] (see 11). 

The quantity needed to deal with time-dependent phenomena in DFT is 
the quantal action M. According to WKG, this is 



Jto 2 J 



dr 



^t(r,t)^^ a^^^^'^) 



-H> . (2) 
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Following the well-known DFT argument, we shall prove below that the 
potentials are in one-to-one correspondence with appropriate gauge-invariant 
densities (Theorem I); that the action functional can be written in terms of 
these densities (Theorem II); and that a practical scheme can be given to 
map the interacting system into a non-interacting one driven by effective 
potentials which include exchange and correlation (Theorem III). 

The action (|]) is invariant under the gauge transformation A(r, t) — »• 
A(r,t) + hVA{r,t) and V{r,t) V{Y,t) - hdA{r,t)/dt, where A(r,t) is a 
real scalar function such that A(r, to) = mod (27r). The source function 
rj{r,t) transforms into r]{r,t) exp[iA(r, t)], so that S*^ is gauge- invariant. The 
operators transform according to 

^(r,t) ^^(r,t)exp[^A(r,t)] (3) 

and 

j(r, t) ^ j(r, t) + ih/m)^^r, t)V^(r, t) VA(r, t) , (4) 

with j the paramagnetic current- density operator. We choose the following 
gauge-invariant densities: the current density j(r, t) = (j(r, t))—n{r, t)A(r, t)/m, 
the condensate density nc(r,t) = \{ip'>{r,t))\'^ and the velocity field Vs(r,t) = 
{h/'m)'Vip{r,t) — A(r, t)/m. The total density n{r,t) is not an independent 
quantity, since it is determined by the continuity equation. The pair of physi- 
cal quantities nc(r, t) and Vs(r, t) can be replaced by the gauge- invariant order 
parameter $(r, t) = {ip{r,t)) exTp[—{i/h) fl:^dt'V{r,t')] : we shall exploit this 
fact whenever convenient. Finally, we stress that is the irrotational part 
of the velocity field. 

The role of A and V . We have already discussed the role of the source 
field ?7, which is characteristic of a neutral superfluid. Before proceeding 
to prove the TD-DFT theorems we pause to discuss the physical import of 
the fields A and V entering the Hamiltonian (|l]). Their meaning is obvious 
for a superconductor, but needs elaboration for a superfluid. We use as a 
guideline for this purpose the linearized two-fluid model below threshold for 
vortex formation. 
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We consider first the vector potential A. As pointed out by Baym while 



discussing the rotating-bucket experiment [|T5|, what may create transverse 
currents in a superfluid is a spoon-stirring mechanism. Let uj be the stirring 
angular velocity, with magnitude below the threshold lOc for vortex formation. 
If L = mr X j is the angular momentum, the corresponding Hamiltonian term 
can be written as / c/r ■ L = / (ir j ■ (a; x r). Comparing with the minimal 
coupling form j ■ A in the Hamiltonian the component of A parallel to 
j is A = mo; x r, namely with m times the rigid-body rotational velocity 
of the fluid. Since that part of the fluid which can respond to a transverse 
probe is by definition the normal- fluid component, lj x r is the normal-fluid 
velocity v„ and thus A = mv„. This result remains true for a non-rotating 
fluid, as demonstrated by Hohenberg and Martin by means of a Galileian 



transformation 14 



Let us turn to the scalar potential V . The quantity V + hd^p/dt^ with 
(/9 being the phase of the condensate, is gauge-invariant. Therefore, writing 
the equations in the gauge in which the scalar potential vanishes corresponds 
to a Galileian transformation to a reference frame moving with velocity v^. 
This fact will be used in the proof of Theorem I below. 

From the above arguments regarding the potentials A and V , it follows 
that d{wn — \'g)/dt = dA/dt + VV^/m and therefore is gauge- invariant. In 
fact, in the two-fluid model (with n = ps + Pn, Ps and p„ being the super- 
and normal-fluid densities) the gauge-invariant current density is = p^v^ -|- 
PnVn — nvn = Ps{^s — v„). This is the current as seen in a reference frame 
which moves with the normal-fluid component and determining one of the 



driving forces in the Landau-Khalatnikov equations . 

We conclude by remarking that a parallel can be made between the two- 
fluid equations for neutral superfluids and Maxwell's equations for charged 
superconductors. From the above analysis it turns out that the equation 
Vx(E+c"^9A/(9t) = or else E+c~^ dA/dt = —VF is just the condition for 
irrotational flow. The "electric field" E is identified with {m/ ps)djr/dt. As 
expected, the Maxwell equation for V x B expresses the relation of continuity 
between particle and current densities. 
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After this excursus we return to the basic theorems of TD-DFT for neutral 
superfluids. 

Theorem I. It states that the densities {d} = {j(r, t), $(r, t)} are uniquely 
related to the potentials {p} = {A{r,t),V{T,t),ri{r,t)}. One has to show 
that two sets of potentials {p} and {p'}, which differ by more than a gauge 
transformation and can be expanded in Taylor series around to, determine 
two different sets of densities {d} and {d'} evolving from a common initial 
equilibrium state. While the statement is trivially true at time Iq, it is 
sufficient to prove it at some time t infinitesimally later than to by relating 
the coefficients of the Taylor series for the densities to those for the potentials 

i- 

We thus consider the Heisenberg equations of motion for the densities. 
The potentials contribute to the equation for the induced current ja with 
terms including nVaV, ApV pja-, j/sVaAjS, j^V pAp^ and {hn/m)Ai3V aAp. It 
is evident that the proof will be easier in reference frame moving with velocity 
v^: in this gauge, as already remarked, the scalar potential vanishes and j 
and A are both transverse, so that all the above terms vanish. We proceed 
within this gauge, signalled henceforth by a tilde on the potentials. 

Since the two sets of potentials are different, their Taylor-expansion co- 
efficients must differ at some order, say I for A and A' and /' for fj and fj' . 
It is then sufficient to show, for the lower among / and Z', that different co- 
efficients in the expansion of the potentials imply different coefficients in the 
expansion of the densities |§-|[T^. In the case I < I' we have 



|r [j(r,t) - j'(r,t)],.,„ = [A(r,t) - A'(r,t)]^^^^ , (5) 

while in the case I > I' we have 

'^1^ ["^(r, t) - $ '(r, = 1^ [^(r, t) - r^'(r, t)],^,^ . (6) 

As a consequence of eqs. and (^, the set of densities {d} will differ from 
{d'} at times infinitesimally later than to. Hence they are different. This 
proves Theorem I. 
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The conclusion thus is that in a superfluid the potentials are unique func- 
tionals of the densities. Since from the Heisenberg equations of motion the 
field operators are functionals of the potentials, we may state that the ensem- 
ble expectation value of any gauge-invariant operator is a unique functional 
of the chosen set of densities. 

Theorem II. It states that (i) the action Q in given external potentials can be 
expressed as a unique functional [{d}] of the densities {d} where the super- 
script indicates the external potentials; and (ii) [{d}] is stationary with 
respect to the actual densities {d^} of the interacting system. The proof pre- 
cisely parallels that given by WKG for superconductors, once their complex 
gap function A(r, t) is replaced by $(r, t) or by the subset {nc(r, i), Vs(r, t)}. 
The functional is given by 



Q° [{d}\ = R [{d}\ - W [{d}\ - P° [{d}\ - [{d}] 



(7) 



where 



R [{d}] = (1/2) /* dt' I dr{jjj^ [{d}] \{ihd/dt') - (ny2m)V^] ^jj [{d}]) + c.c. 

Jto J ^ 

(8) 

and W [{d}] = dt'{W [{d}] (f)) are its universal parts, while 



[{d}] ^ £ dt' I dv [(V%r, t') + ^A'\r, t')) n [{d{r, t')}]+ 
+ AO(r, t') ■ (j(r, t') - n [{d(v, t')}]A [{d(r, t')}] /m) 



and 



to 



S'[{d}]^ / dt' dr\^^(r,t')^*{r,t')+vTir,t')^r,t'] 



(9) 



(10) 



depend on the external potentials. The gauge has been chosen so that the 
functional V [{d}] equals the external scalar potential V^{r,t) and r]g is de- 
fined by 7]^g = 7]^exp[-{i/h) dt'V^iv^t')]. 

Along with the basic idea underlying DFT, Theorem II admits a map of 
the densities in the real system onto those of a reference system subject to 
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appropriate potentials. This map is proven in Theorem III, which defines 
the so-called Kohn-Sham scheme needed to implement TD-DFT. 

Theorem III. It states that there exist unique reference-potential functionals 
{p^ [{<^}]} such that the densities {d^} calculated within the chosen reference 
system coincide with the densities {d^} of the real interacting system. 

Following again WKG, we first define the action functional [{d}] = 
[{d}] — [{d}] — [{d}] for the reference system as in eq. (0) for its 
interacting analogue. The functional [{d}] is written as 



Q" M] = R"" [{d}] - [{d}] - [{d}] - [{d}] 



thereby defining the exchange-correlation functional Q^c [{d}]- We can now 
exploite the stationarity of both [{d}] and [{d}] to obtain a set of 
equations relating the potentials {p^ [{d}]} to the original external potentials 

{/}■ 

The resulting equations, in addition to [{(i°(r, t)}] = V^{r,t) are as 
follows: 



(5P^ [{d}] 



SP' [{d}] 
5j(r,t) 



+ 



J{dO} 



5Q,c [{d}] 
5j(r,t) 



(12) 



J{dO} 



and 



<(r,t) + 



SP"" [{d}] 



{dO} 



r,t) + 



SP'[{d}] 
(5$*(r,t) 



+ 



5Q,c [{d}] 



(13) 

with its complex conjugate. These equations define the effective exchange- 
correlation potentials. 

We conclude this discussion by noticing that eqs. (|12D and (|13D have 
been derived in a previous paper []T^ within a linear-response formulation 
of TD-DFT for superfluids. In brief, by writing the microscopic equation 
of motion for the order parameter in terms of the condensate self-energy we 
proved that the matrix expressing the linear response of nc{r,t) and Vs(r, t) 
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to the symmetry-breaking field rj explicitly has the Hohenberg-Kohn-Sham 
structure. 

Reference system. A suitable choice of the reference system for a superfiuid 
at finite temperature is provided by the gapless Bogolubov-Popov approxi- 
mation. This accounts for the thermally excited non-condensate cloud and 
satisfies the Hugenholtz-Pines theorem [|1^. In this approximation the den- 
sities can be written as 

n,(r,t) = |$(r,t)|2, (14) 
v,(r,t) = {h/m)Vip{r,t) - A^(r,t)/m (15) 

and 

j(r, t) = jc(r, t) + j(r, t) - n(r, t)A^(r, t)/m . (16) 
Here jc(r, t) = ?2c(r, t)vs(r, t) is the condensate current and 

j(r, = ^ E [NnUnir, t)VU:{r, t) + (iV„ + 1) K(r, t)VV:{r, t) - c.c] 

(17) 

is the current carried by the non-condensate. In eq. (|T^ f/„ and Vn are the 
Bogolubov functions and Nn = [exp{En/kBT) — is the boson thermal 
factor, with En being the energy eigenvalues in the Bogolubov-Popov equa- 
tions (see below) at the initial time to. Finally, ^(r,^) = nc{r,t) +n{r,t), 
with 

n(r,t) = E K(|f/n(r,t)P + |K(r,t)|') + |K(r,t)|'] (18) 

n 

being the non-condensate density. The anomalous density is given by (V^(r, t)-i^{r, t)) 
$'(r, t) + En (2iV„ + 1) f/„(r, t)V:ir, t). 

In eq. (^) the condensate wave function $(r, t) satisfies the Schrodinger 
equation p 



^^^^|^=/:^$(r,t)+<(r,t), (19) 
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-(2m)-i [^V - iA^(r, t)] + V^{r, t) and 1/^(r, t) = V%r, t) + 
r')n{r',t) + Vxc{r,t), with Vxc{r,t) being determined from eq. 



where : 
2 / dr'w{r 

(P^ . The gauge-invariant reference source field is rj^{v,t) = rpg{ 
$(r,t)/dr'w(r-r')|$(r',t)|2 + r/,,(r,t) = 770(r,t) + 5Q,,/<5$* (see eq. {^). 
In the special case of a point-contact interaction eq. (|19]) becomes the well- 
known Gross-Pitaevskii equation. 

The Bogolubov functions Un and Vn satisfy the single-particle coupled 
equations 




/c/r'w(r-r')$^(r',t) 



/ dY'w{ 



*2/ 



r',t) 



(20) 

We point out that, as a result of imposing gauge invariance, the reference 
system in eq. ([T9|) is the same as that in eq. (|20|) . 

In summary, we have demonstrated the basic Hohenberg-Kohn-Sham- 
type theorems underlying TD-DFT for inhomogeneous neutral superfluids 
at finite temperature below threshold for vortex formation and proposed 
an implementation based on a reference system described by the Bogolubov- 
Popov theory. We have also explicitly pointed out similarities and differences 
with respect to charged superconductors as treated by Wacker et al. [llOi and 



briefly remarked on the linear-response limit as treated by Chiofalo et al. [|T^ 



A final comment is in order. For super-critical rotational velocities quantized 
vortices will appear in the superfluid: at that point the velocity field 
ceases to be irrotational and acquires a regular contribution describing 
the velocity of each point of a vortex line as well as a singular contribution 
due to the quantized structure of the vortex line |TP|. The regular term leads 
to the well-known Magnus force on a vortex line [pOl and to friction forces 



between superfluid and normal-fluid components, which are proportional to 
Vr — v„ [|I9| . Therefore, in order to account for vortices the present TD-DFT 
approach will need extension to include one additional external field and one 
additional density variable. We hope to return to this problem in the near 
future. 



Un{j,t) 

K(r,t) 
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